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Abstract 



We derive the kinetic equations for both the covariant and equal-time 
Wigner functions of Dirac particles with electromagnetic, scalar and pseu- 
doscalar interactions. We emphasize the constraint equations for the spinor 

components in the equal-time formulation. 
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Typeset using REVT[t;X 



One of the recent progresses in transport theory [||,|| is the estabhshment of transport 
equations for spinor ||^ and scalar equal-time Wigner functions with abelian gauge in- 
teraction. The main advantage of the equal-time formulation lies in the fact that the initial 
value of the equal-time Wigner function can be directly obtained from the corresponding 
initial field operators, since there is only one time scale in the equal-time formulation. There- 
fore some quantum problems such as pair production in strong electric fields have so 
far only be solved in the equal-time formulation P,^. The advantages of the covariant 
formulation, on the other hand, are explicit Lorentz invariance and the feature that the 
kinetic equations naturally split into a transport equation of Vlasov-Boltzmann type and a 
generalized mass-shell condition which makes explicit the off-shell effects generated by 
the collisions in the system. In Ref. P] we discussed the relationship between the covariant 
and equal-time kinetic equations. By taking the energy average of the covariant equations 
we obtained both transport and constraint equations for the equal-time Wigner functions in 
scalar and spinor electrodynamics. 

Recently, Shin and Rafelski |^ discussed a system which, in addition to the electromag- 
netic vector interaction, included also scalar and pseudoscalar (strong) interactions. For the 
system defined by the Lagrangian 



they derived a generalized group of transport equations for the spinor components of the 
equal-time Wigner function in the mean field approximation. Here A^, a and vr are elec- 
tromagnetic, scalar and pseudoscalar fields, respectively, g^j and g^, are the scalar and pseu- 
doscalar coupling constants, and UMicr,'^') is the mesonic self-interaction potential. In this 
short paper, we first study the full covariant kinetic equations for such a system and then 
derive from them transport and constraint equations also for the equal-time Wigner func- 
tion, using our recently introduced energy averaging method. We thereby supplement the 
equal-time transport equations derived in Ref. by a set of non-trivial constraint equa- 



C = ^p{i-f^'^^ - e^^A^ -m + g„a + ig^'K-i^)^) - ^F^^F^" 
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tions. The importance of the latter is emphasized both for the semiclassical hmit and for 
the general quantum dynamics. 

In the following we restrict ourselves to the kinetic equations for the spinor Wigner 
function in the mean field approximation, by replacing the field operators A^, a and tt by 
their mean values. Such a Hartree approximation can be justified for sufficiently strong 
fields and has so far been used in most applications of quantum transport theory. Since 
the meson potential Um does not explicitly appear in the kinetic equations for the spinor 
Wigner function (it only enters the selfconsistency equations for the scalar and pseudoscalar 
mean fields), the discussions below holds universally for all potentials WAf(o", tt). 

The covariant spinor Wigner function Wi{x,p) is the ensemble average of the Wigner 
operator W^^x^p) which is the 4-dimensional Wigner transformation of the covariant density 
operator <i'4(x, y): 



Calculating the first-order derivatives of the covariant density matrix $4(0;,?/) with re- 
spect to X and y and using the Dirac equation 



one obtains an evolution equation for $4(0;,?/). After Fourier transforming with respect to 
the 4 components of the relative coordinate y one obtains for arbitrary external fields A^, a 
and TT the generalized VGE [|lOl equation 




(2) 



^lii\px + «ey4^(x) = m'ijj{x) — go-(T{x)'ijj{x) — ig.„TT{x)'y^'ijj{x) 



(3) 



(4) 



where now 



M 



Ml + iMr 



2 , 
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^t^ix.p) =p^-ieh / ^dss Ff,^{x - ihsdp) 



Df,{x, p) = d^ - e Ids F^y{x - ihsdp) . 



D^{x,p) = g^cos ^(^) ' 

Ml {x, p) = m — g„ cos ^2^^ ^^""^^ 

M2{x,p) = gasin ^ 
A = d^-dp . 



(5) 



For the discussion of the semiclassical expansion below we have exphcitly given all factors 
of h. If one neglects in these equations the abelian gauge field tensor and all terms of second 



or higher order in h one recovers the results recently obtained by Florkowski et al. [ITT 



Inserting the spinor decomposition |T0 

W4ix,p) = 



F{x,p) + ii^P{x,p) + 7^\/^(x,p) + 7^75A''(x,p) + -(J^^S^'''{x,p) 



(6) 



we split the complex equation (H) into 10 real equations for the spinor components: 



WV^-D^P = M^F, 
hD''A^ + 2D5F = 2MiP, 
U^F + ^D''S,,-U,A^ = M,V, 



-fiD^P + e^.^.WS'^f - 2U,V^ = 2MiA, 

hD^V^ + 2Y{^P = 2M2F , 
WA^ + U,F = -M2P, 
^D^F - WS,^ - D,A^ = M2V; , 

ft 1 



(7) 



Again the corresponding results up to first order in h can be found in Ref. |TT[ if we do not 
consider the abehan gauge interaction. 

Following the treatment of Ref. [Q, we write the equal-time Wigner function Wslx, p) as 
the energy average of the covariant Wigner function W4{x,p) and decompose it in a similar 
way [0 into its spinor components: 



Wsix, p) = J d^ye-'^ y (^x + |, exp te J_\ rfsA(x + sy, t)-y 

dpoW4{x,p)'yo 



= I [fo{x, p) + 75/1(3;, p) - 27075/2(2;, p) + 70/3(3;, p) 

+757o7-go(a;, p) + 7o7-gi(a;, p) - il-g2ix, p) - 757-g3(a;, p) 

The second line in (H) suggests that the equal-time kinetic equations may be obtained by 
taking the energy average of the corresponding covariant equations. As pointed out in Ref. 
P] and seen below, this procedure yields additional information on the equal-time Wigner 
function, in the form of additional constraint equations which are lost in the procedure used 
in Ref. which is based on a 3-dimensional Wigner transformation of the equations of 
motion for the equal-time density matrix. 

The energy average of the equations (|^ gives rise to kinetic equations for the spinor 
components /» and gj (i = 0, 1, 2, 3) of the equal-time Wigner function. In direct analogy to 
Ref. one obtains a set of transport equations which are identical to the ones derived in 
i, 

^(A/o + D-gi) = -2g„(Tof3 + 2g^nof2 , 
^(A/i + D-go) = 2g^(7ef2 + 2gn'n-ef3-'2mf2, 
^(A/2 + 2n-g3) = ~2g^crJi + 2g^7rofo + 2mfi, 
^(A/3-2n-g2) = -2g^aofo-2gnTrefu 
h{Dtgo + D/i) - 2nxgi = -2g^(Tog3 + 2g^TTog2 , 
^(Agi + D/o) - 211 X go = 2g^aeg2 + 2g^TTeg3 - 2mg2 , 
^(Ag2 + Dxgg) + 2n/3 = -2g^aegi + 2g^TTogo + 2mgi, 



h{Dtg3 - D X g2) - 2n/2 = -2g„aogo-2g^7regi , 



(9) 



as well set of additional constraint equations: 



dpopoF - -hB ■ g2 + flo/s - gn-^ofi = (m - gaCre)fo , 

dpoPoP + ■ g3 + 60/2 - S'TrTTe/o " fl-aO-o/l = , 

dpoPoVo - n ■ gi + fio/o - 9nTTef2 = (m - g^ae)h , 

dpoPoV - ^/iD X go - n/o + flogl + S'vrTl'ogS + fi'aO-og2 = , 
dpoPoAo + n ■ go - fio/l - S'ttTTo/s - 5'aCro/2 = , 

dpoPoA + i/lD X gi + n/i - flogo + fi'7r7reg2 = -(m - fi'<70-e)g3 , 
dpoPoS^'Gi - + n X g3 - nog2 + S'TrTTegO + S'aCTogl = , 

j dpopoSjke^'''ei - h'Df2 + 211 x ga + 2nog3 - 2g^TTogi = 2(m - gaCre)go 
Here we defined the following 3- dimensional operators 



Dt 
D 

n 

no 

(To 
TTe 



x,p 

X, p 

x,p 

x,p 
x,p 
x,p 
x,p 



dt + e ^ rfsE(x + zsWp,t)-Vp, 
V + e rfs B(x + ishVp, t) x 
p — iefi Idss B(x + ishVp, t) x Vp . 
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h / ^ (issE(x + isWp,t)-Vp. 



cos(-V ■ Vp)a(x) , 
sin(^V ■ Vp)(T(x) , 
cos(^V ■ Vp)7r(x) , 
sin(^V ■ Vp)7r(x) , 



(10) 



which correspond to the covariant operators (^. E and B are the electric and magnetic 
fields, respectively. 

We now discuss the relevance of the additional constraint equations (pUD, both in the 
classical limit and for genuine quantum situations. In the classical limit, the covariant kinetic 



equations have solutions P] of the form W^{x,p) = W^{x,p) 5{po =F Ep). Here Ep is still 
an arbitrary function of p; its concrete form will be seen momentarily as a result of the 
constraint equations. To zeroth order in h, the kinetic equations (|]) and (0) lead to the 
following constraints for the classical spinor components and the quasiparticle energy Ep-. 



Ep 
Ep 



ft = 



Ep 



Ep 



P X + Kg^ 
m — Va 

^ E^jm- K)g^ - (m - VMp-So) T EpV^p x 

Ep{m*y 



E'^ = {m*y + p'. (12) 
Here we defined 

V. = g.a, V^ = g^7i, (m*)^ = (m - K)' + Kr' • (13) 

The classical transport equations for the charge and spin densities are seen to originate from 
the equation (|^) at first order in h: 

dtft + d gf = iF^-Vp)ft - (F.-Vp)/2± , 







where 

dt = dt + eE- Vp , d = V + eB X Vp , 

= -f7.Va , = -g^VTi . (15) 
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Introducing the particle and antiparticle charge densities / and / and spin densities g and 
g through 

f{x, p) = fo{x, p) , /(x, p) = /o"(x, -p) , 

g(x,p) = g([(a;,p) , g(x,p) = g(^(a;, -p) , (16) 

and using the constraints (jl^j), one derives from ( ]T^ ) a Vlasov-type transport equation for 
/, 

dj + (v-V)/ + (^eE + ev X B - j .^j ^ , (17) 

and a similar equation for the spin density g: 

^tg + (v- V)g + I eE + ev X B 1 ■ Vpg 

= ^((p.g)E-(E.p)g)-^Bxg 

, (m - K)F^ + KF. / (p-g) \ (m - K)F, - KF^ , . 

+ — — 'i^^^^j^ — ^s^^^^? — '^^'^^ 

, {m-K)dtV^ + V^d,V / (p-g) K \ 

+ (m - K)^(m*)2 ~ - ^ 

((P-g)P + (- ) gj + (;;^P X g - ^gj . (18) 

Here v = is the velocity of the classical transport flow. The corresponding equations for 
the antiparticle densities / and g are obtained by changing the sign of the electric charge, 
e — i> — e. The spin evolution equation (p!8| ) is the 3-dimensional phase-space version of a gen- 
eralized Bargmann-Michel-Telegdi (BMT) equation ]T^pD| which describes spin precession 



in external electromagnetic, scalar and pseudoscalar flelds. The flrst two terms of the r.h.s. 
which couple the spin density to the electromagnetic fleld have already been given in 
The remaining terms describe spin-precession effects caused by the coupling to the scalar 
and pseudoscalar flelds. They are the 3-dimensional version of the much more compact 
covariant expressions derived in [jll . 
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It is well known [|TO|JTT| that the covariant spin equation can be greatly simplified by 



introducing the covariant "spin up" and "spin down" distribution functions 

F^,{x,p) = ^± ^^^'^P^f'^'^P\ (19) 
m — Vo- m*[x) 

where S*^ = A^/{—A-AY^'^ is the covariant and normalized spin phase-space density. 
The kinetic equations for the analogous equal-time distribution functions f±s{x,p) = 
EpJdpoF±s{x,p)5{po - Ep) and /±,(x, p) = Ep J dpo F±six, -p) S{po + Ep) decouple into 
scalar equations of the form (|1^ also obeyed by the charge density: 

dtf±s + (v-V)/±. + (^eE + ev X B - " " ^^^. ^ ^ ^ ^ ^^q) 

and a similar one for the antiparticle distribution f±s obtained through e — e. 

Although the energy averaging method and the procedure of Ref. lead to the same 
Vlasov transport equations for / and /, we would like to point out an essential difference: 
Since the equal-time constraint equations (|10D do not disappear even in the classical limit, 
they reduce (together with the transport equations (^) to order h^) the number of indepen- 
dent spinor components from 8 in ||^ (namely /o, /s, go and ga) to 4 here (namely the charge 
density /o and the spin density go). Furthermore, the form of the quasiparticle dispersion 
relation Ep need not be guessed as in 0, but follows from the constraints. In the treatment 
of Ref. 0, the equations for the charge density /o and the mass density /s remained coupled, 
and the Vlasov-type transport equation (|T7|) was only obtained after a suitable redefinition 
of the classical distribution functions by hand (see equations (57) and (58) in [§). Our 
constraint equations (0) turn these redefinition equations into identities. 

For genuine quantum problems like pair production, one must consider the non- 
perturpative (in Ti) form of the constraint equations (0). To get a closed set of equations 
for the spinor components of the equal-time Wigner function, one eliminates the first p^- 
moments / dpopo W4^{x,p) from ([T0|) by combining them with the first po-nioments of those 
covariant equations in (^ which gave rise to the equal-time transport equations (H). The 
procedure is entirely analogous to the one detailed in Ref. ^ , and one obtains the constraints 

L/o + M-gi-F.„-g2-F^„-g3 = 0, 

9 



L/i + M ■ go - (F,, + F,) • g2 + (F^, + F,) ■ g3 = , 
Lf2 + 2N ■ g3 - (F^, - F^) ■ gi - F^„ • go = , 

Lfs - 2N • g2 + (F., - F.) • gi - F^„ ■ go = , (21) 
Lgo - M/i - 2N X gi - F,j2 - F,j3 - (F., - F^) x ga - (F,^ - F,) x g2 = , 
Lgi - M/o - 2N X go - (F^^ - F^)/2 + (F.^ - F.)/3 + F,„ x g3 - F^„ x g2 = , 
Lg2 + M X g3 - 2N/3 + F.Jo + F^Ji - (F,,^ - F.) x go + F^„ x gi = , 
Lg3 - M X g2 + 2N/2 + F, Jo - (F,, - F.)/i - (F,, - F,) x go - F.„ x gi = , 

with 

1 

L{x, p) = ie J^^ ds s (y X B(x + isVp, t)) ■ Vp , 



2 

. 1 



X, 



M(x,p) =ie f\ rfss v(E(x + isVp,t)-Vp) +e f\ ds (E(x+zsVp,t) - E( 
-'-1 ■'-1 

1 

N(x, p) = -e Ids (Vp-V) E(x + isVp, t) 

+e / ds s^{dtB{x + isVp, t)) xVp + ie ^ (issE(x + isVp, t) , 

Fa„ = -ga^(To , 
Fff, = -Qa^de , 
FtTo = -S'tt VTTo , 

F., = -^.Vvr, . (22) 

These constraints, which are closely associated with the quantum corrections to the clas- 
sical mass-shell condition, hold for arbitrary external fields. Only for homogeneous fields, 
the spatial derivatives of the electric, magnetic, scalar and pseudoscalar fields vanish, and 
therefore the constraint equations (^Tj) disappear identically. 

These results extend our work in |^ where we studied only electromagnetic interactions. 



They also complement the recent studies of Refs. |Tl[ and by generalizing the former 
to all orders in h and supplementing the latter with the necessary additional constraints. 
The generalization to scalar and pseudoscalar interactions is relevant for the development 
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of a chiral kinetic theory with possible apphcations to the non-equihbrium dynamics in 
relativistic heavy ion colhsions. Eq. (0) has the generic form of an effective chiral Lagrangian 
for low-energy hadronic physics; suitable choices for the meson potential Um correspond to 
different popular effective chiral models. For example, replacing the second line in Eq. (0) 
by C{a^ + vr^) [|rT| yields the semi-bosonized version of the Nambu-Jona-Lasinio model 
[[13| . On the other hand, the meson potential Um ~ (o"^ + tt^ — f ^)^, with t> as a parameter, 
represents the original linear a model Jill ■ These two models are widely used to describe low 



energy hadronic physics and the chiral phase transition at high temperatures and densities. 
Therefore, the study of the kinetic equations in both covariant and equal-time formulations 
for the Lagrangian (|lD is helpful for an understanding of the dynamical consequences of the 
chiral properties of QCD. 

We should not close without repeating from Ref. [§] an important cautionary remark. 
The energy averaging procedure elucidates that the kinetic equations for the equal-time 
Wigner function W2,{x, p) = / cipo Wis^{x,p) form only the lowest level of a coupled hierarchy 
of equations for the energy moments of the covariant Wigner function, / dp^p^ Wi^{x,p), n = 
0, 1, 2, . . . The full hierarchy can be obtained by taking the energy average of the appropriate 
Po" moments of the generalized VGE equation (^). Its properties are presently under study 
in the context of the simpler case of scalar QED ||15|. The connection between certain 



truncation schemes for this hierarchy of equal-time moment equations and the well-known 
gradient expansion in the covariant approach |l|,0,|l^ is an interesting subject for further 
research. 
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